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ABSTRACT 

We present the first large-scale stellar kinematic and line-strength maps for dwarf elliptical 
galaxies (9 in the Virgo Cluster and 3 in the field environment) obtained with the SAURON 
integral-field unit. No two galaxies in our sample are alike: we see that the level of rotation 
is not tied to flattening (we have, e.g. round rotators and flattened nonrotators); we observe 
kinematic twists in 1 Virgo and 1 field object; we discover large-scale kinematically-decoupled 
components in 2 field galaxies; we see varying gradients in line-strength maps, from nearly 
flat to strongly peaked in the center. The great variety of morphological, kinematic, and stellar 
population parameters seen in our data points to a formation scenario in which properties are 
shaped stochastically. A combined effect of ram-pressure stripping and galaxy harassment is 
the most probable explanation. We show the need for a comprehensive analysis of kinematic, 
dynamical, and stellar population properties which will enable us to place dwarf ellipticals and 
processes that govern their evolution in the wider context of galaxy formation. 

Key words: galaxies: dwarf - galaxies: evolution - galaxies: formation - galaxies: kinematics 
and dynamics - galaxies: populations/stellar content 



1 INTRODUCTION 

The importance of dwarf elliptical (dE) galaxies lies in numbers. 
They are the numerically dominant galaxy class in clusters, by far 
outnumbering any other type (Sandage et al. 1985, Binggeli et al. 
19881. They are low-luminosity objects (Nig ^ — 18 mag), which 
makes both external influences and internal feedback mechanisms 
far more extreme than in more massive galaxies. Thus, they are 
ideal "laboratories" to test various mechanisms that shape galax- 
ies. They are also a surprisingly (from a historical point of view) 
inhomogeneous class: which has made it challenging both to relate 
different dE subtypes to each other, as well as to place the whole 
class in the larger context of galaxy assembly and (trans)formation 
processes. This is perhaps best reflected by the fact that even today 
the astronomical community is not unanimous in their choice of the 
most accurate name for this galaxy class. 



1.1 Photometric, stellar population, and kinematic properties 

Much has been learned about dEs in recent years. As far as pho- 
tometric studies go, the first indications for a disky nature of some 
dwarf early-type galaxies are found in Sandage & Binggeli ( 1984) 
and Binggeli & Came ron] ( |1991^ . The first pieces of evidence for 
the presence of substructures such as spiral arms were found by 
|Jerjen et al"] ( |2000fr and |Barazza et H\ ( |2002| > for Virgo, |Graham| 
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eTaT1p003) for Coma, and |De Rijcke et al.| < |2"003l > for Fornax dEs. 
Lisker et al. (20061 were the first to carry out a systematic study 
of all Virgo dEs down to mg=18 mag covered by the Sloan Digital 
Sky Survey Data Release 4. Their work resulted in a new classifica- 
tion that divided dEs into 4 sub-types: nucleated (N), non-nucleated 
(nN), blue-core (be), and disky (di) objects. They obtained the nu- 
merical contribution of each subtype to the whole class, discovered 
that the level of substructure scales with galaxy mass, as well as 
that nucleated galaxies are preferentially found in the cluster's cen- 
ter. Their results also suggested that a dichotomy exists in the dE 
group that is based on the presence of a nucleus. However, the ACS 
Virgo Cluster Survey (Cote et al. 2006 1 found nuclei where they had 
not been seen previously, suggesting that the presumed dichotomy 
might have arisen from data quality limitations and that there might, 
in fact, exist a continuity in the nucleus-to-galaxy body ratio. With 
their deep photometry of 100 Virgo dwarfs, [Tanz et al.| (2012) estab- 
lished that the level of substructure was previously greatly underes- 
timated (only 25% of their galaxies are fitted with a single Sersic 
function). 

Similar diversity of properties of dEs has been found through 
the analysis of their stellar populations. dEs display a variety of ages 
that seem to be uncorrelated with the luminosities of their main bod- 
ies, though nuclei are on average younger and more metal-rich than 
the main bodies (Paudel 201 1) and their metallicity (Z) does corre- 
late with luminosity. There also seems to exist a relation between 
the ages and the projected distances of dEs to the Virgo Cluster's 
center ( Michielsen et al. 2008). As for the gradients in these prop- 
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erties, |Chilingarian| p009) and |Koleva eTal] p009) find that dEs 
show a large variety of metallicity gradients; the latter authors also 
note the presence of flat metallicity profiles in disky galaxies. Kol- 
|eva et al.] ( |2011| > later find no relation of Z gradients with luminos- 
ity, velocity dispersion, central age, or age gradient, but they do find 
a link between age gradients and luminosity, with fainter galaxies 
(Mb >-17 mag) having stronger positive gradients. 

Kinematic studies confirm and add to the variety. Rotating dEs 
are preferentially found in the cluster's outskirts and non-rotating 
in its center fTbloba et al.|2009F . |Toioba et al.| ( [20TT] > find that ro- 
tationally supported dEs have rotation curves that resemble those 
of late-type spirals. The studies of |Geha et al.| (2002) and|Toloba 



et al. 



2006 



( |201 1| > also suggest that dEs (like the Es of Cappellari et al. 



but unlike the dSphs of Wolf et al. 2010} are not dark-matter 



dominated within one effective radius. 



1.2 Relation to other galaxy classes and formation scenarios 

Determining the relation of dEs to other galaxy classes has long 
been a difficult task. Their seemingly simple and round appearance 
caused them to be classified together with giant ellipticals (hence 
the name) as their low-mass counterparts. A general consensus is 
yet to be reached regarding the (dis)continuity and the physical 
meaning of various scaling relations that include the dE, E, and SO 
classes (e.g. |Graham & Guz man 2003 Ferrarese et al. 2006| [Kor^| 
|mendy et al.|2009) . 

However, once the variety of properties of dEs became clear, 
the scenario in which dEs are primordial objects which expelled 
their gas in early stages of their evolution because of supernova 
explosions (e.g. |Mori et al7|| 1999[ > could only hold for the round, 
pressure supported systems found in the central parts of the cluster, 
since they are the only ones that seem to have only old populations, 
with no star-formation episodes following their initial creation. The 
varying range of properties discussed above strongly suggests that 
the majority, if not all, dEs are most likely transformed late-type 
(spiral and irregular) galaxies that entered the cluster at later stages 
and had their properties shaped by external perturbations. This 
notion has, in fact, been around for almost three decades. The 
first authors to suggest it were |Wirth &~ Gallagher (1984) and 
|Kormendy] 1 1985|>, and it wa s later reiterated by Binggeli et al.| 
T988) , |Binggeli & Cameron|(T99T) , and |Bender et al.|(1992| >. 



Different scenarios have been proposed in which star-forming 
late-type galaxies could have been transformed into the quiescent 
dEs we see today, the main ones being: 

(i) ram-pressure stripping, i.e. interactions between galaxies 
and the intergalactic medium (Gunn & Gott|[T972"] |Lin & Faber| 
1983), transforms gas-rich dwarf irregulars/spirals into galaxies 
with very little cold gas, thus effectively stopping star formation. 
In this case it is expected that the intrinsic angular momentum of 
the progenitor should be conserved. 

(ii) galaxy harassment, i.e. multiple high-speed encounters 
with other galaxies in the cluster environment, combined with tidal 
heating resulting from the interactions with the cluster's potential 
well (Moore et al.|1998|[Mastropietro et al.|2 005 1, can remove stel- 
lar mass and significantly change the morphology, the strength of 
the effect likely being tied to the local environment density. In this 
case, a galaxy is expected to lose some (but not all) of its intrinsic 
angular momentum. 

There exist examples of environmental transformation pro- 
cesses in action. We know that galaxies are falling into the Virgo 
Cluster and are having their gas stripped from, for example, the 
results of Chung et al. (2007) who show spiral galaxies with Hi 
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Figure 1. Top panel: positions in the cluster for all Virgo dEs from |Lisker| 
|et al. |(2007} with the four (disky, blue core, nucleated, non-nucleated) sub- 
types overplotted, showing their concentration in the cluster. Our sample is 
shown as large filled black squares: it contains 3 dE(di)s and 1 dE(bc), the 
remaining 5 are dE(N)s. M87 is denoted with a filled yellow circle. Bottom 
panel: projected distance from M87 vs. ellipticity for the same sample; the 3 
open black squares at the top accompanied by an arrow are our field galaxies 
included here to show their ellipticities. 



tails pointing away from the cluster's center. The same authors later 
confirm a relation between the size and gas content of spirals as a 
function of clustrocetric distance (Chung et al. 20091. Indications 
for harassment are found in the existence of very low surface 
brightness dwarfs, e.g. VCC1052 of |Sandage & Binggeli| (1984) , 
interpreted to be extreme examples of harassed objects, and in the 
presence of faint, round, non-rotating dEs in the central regions of 
the cluster (see the review in Kormendy & Bender |2012| l. 

The degree to which galaxies are affected by the environment 
depends on the galaxy mass (or, more precisely, the depth of the 
potential well, i.e. the degree to which it is able to gravitationally 
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counterbalance the outside "pull") and the intensity with which is 
being disrupted, i.e. the environment's local density. Therefore, we 
can expect relations between the local density (or approximately - 
clustrocentric distance) and the following quantities: 

• Shape and the degree of substructure. High-speed encoun- 
ters with other cluster galaxies will heat up galaxies, making them 
rounder ( |Moore et al.|1 998 Mastropie tro et al.|2005} , although the 
process is also suggested to be able to produce substructures like 
bars or spiral arms. On the observational side |Lisker et al.| ( |2007} 
show that Virgo dEs containing a substructure (disks, bars, or spiral 
arms) are not centrally concentrated, but the dEs without substruc- 
tures are. We also know from Sabatini et al. (2005 ) and De Rijcke 
et al. (2010) that the late-to-early type dwarf ratio decreases with 
decreasing clustrocentric distance. 

• Degree of ordered over random motion. As gas is (ram- 
pressure) stripped from a galaxy, the system is heated up and so 
the V /a drops. Still, it is harassment that is expected to be the pri- 
mary process responsible for removing large portions of rotation, as 
predicted by Moore et al. ( 1996) and Mastropietro et al. (2005 1 and 
observationally suggested by Toloba et al. l |2009) . 

• Dark to stellar matter ratio. The simulations of IMoore et al.l 
( |1998| > and |Smith et al.| ( |2010| > suggest that (at a fixed radius) dark 
matter might be lost more easily due to its different orbital distribu- 
tion. No observational confirmation of this trend has thus far been 
available. 

• Global stellar population properties. Various environment- 
dependent mechanisms, such as ram-pressure stripping (Boselli 
|et al.|20 08 1 or strangulation ( van den Bosch et al.|2008) , are said to 
be able to quench star formation, though their relative importance 
is not yet known. A number of studies find negative age and/or pos- 
itive metallicity global (i.e. among, not within galaxies) gradients: 
Tortora et al.1([20TT],pKngson et al.| ( |200T] >, |Srriith et al.| < [2009] l, and 



Weinmann et al. (2009) find that environment affects colors/star- 



formation rates; this results from the fact that central galaxies have 
resided longer in the cluster and have thus experienced environmen- 
tal influences for a more extended period of time. 



1.3 This study 

We can see that a substantial amount of effort has recently been put 
into trying to understand the dE class. What is still missing is an 
"all-inclusive", comprehensive analysis of all observable properties 
that is later matched with simulations in order to put qualitative 
constraints on the (trans)formation processes. 

In response to these issues we have started a study of dEs in 
the Virgo Cluster using the SAURON integral field unit (IFU). Us- 
ing integral-field spectroscopy enables us to simultaneously obtain 
spectra at each position in the field of view and study both kinemat- 
ics and stellar populations of the galaxies in two dimensions, giving 
an unprecedented view of the mass distribution and star formation 
histories. The two-dimensional capabilities and large spatial cov- 
erage of SAURON are key to clearly identifying any substructure 
and robustly determining any radial trend in these galaxies. These 
goals cannot be easily achieved with traditional long-slit spectro- 
graphs, where in order to achieve comparable spatial coverage mul- 
tiple pointings would be required. 

This is the first of a series of papers in which we investigate the 
properties of dEs on a sample of 12 galaxies, 9 in the Virgo Cluster 
and 3 in the field. The paper is structured as follows. Our sample se- 
lection, observations, and data reduction are presented in Section 2. 
In Sections 3 and 4 we describe the methods used in the analysis of 



the kinematic and line-strength properties of our sample. Section 5 
presents our results, and the discussion and conclusions are in Sec- 
tion 6. In the other papers of this series we will construct dynamical 
models of dEs, infer their detailed star formation and enrichment 
histories, and compare the global relations of dEs to other galaxy 
classes. All of this together will help us understand and provide 
more insight on the relations listed in the previous subsection. 



2 DATA 

2.1 Sample selection 

The Virgo Cluster is the closest large concentration of dE galaxies 
beyond the Local Group and presents itself as an excellent environ- 
ment in which to study these galaxies. In our target selection we 
used the catalog of |Lisker et al.| ( [2007| ) to select our cluster targets 
and the catalog of the MAGPOP-ITP collaboration (see Mi chielsen] 
|et al.|20 08) to select our field galaxies. 

As for our cluster objects, for our first observing run we 
selected the brightest dwarf galaxies from the Virgo luminosity 
function, attempting to pick objects for which long-slit data showed 
hints of peculiarities. Later, we decided to add objects in such 
a way as to probe galaxies at different values of ellipticity and 
distance from Virgo's central galaxy M87 (see Figure[TJ. This was 
motivated by the assumption that the extent of the environmental 
influence is broadly a function of clustrocentric distance. Thus, 
our objects are in the 0.01-0.56 ellipticity range and 0.78-5.42° 
(0.22-1.56 Mpc) projected distance range. 

The limitations of our attempts at achieving this particular 
parameter coverage are obvious. Projection effects influence both 
the observed flattening as well as the perceived location in the 
cluster. For example, one of our targets,VCC 1261, appears to lie 
close to the center (~ 1.6° in projection), but we later discovered 
from the distance measurements of Mei et al. (2007) that it is, in 
fact, located ~1.6 Mpc behind M87, which translates into their 
intrinsic separation of -~5.8° at the average Virgo distance. 

Thus, any discussion about clustering properties needs to take 
the depth of the cluster into account. For example, objects belong- 
ing to subclasses for which no projected clustering is found (dE(di)s 
of |Lisker et al.|2007| see top panel of Figure [TJ are in fact intrin- 
sically preferentially found in the cluster's outskirts. Also, getting 
individual intrinsic axial ratios requires prior information on the dy- 
namical structure of galaxies (though statistical properties of whole 
subclasses can be assessed, see, e.g. Lisker et al. j 2007). All this 
means that not only the sample size, but also the difficulty in deter- 
mining intrinsic characteristics contribute to the fact that our sample 
cannot be regarded as representative. 

The reason for the inclusion of field galaxies in our sample was 
to enable a comparative analysis of cluster vs. field environments. 
We first preselected objects with which we were able to fill the parts 
of our nights when Virgo could not be observed. From this set we 
picked those with the highest surface brightnesses. We also tried 
to select galaxies for which the closest neighbors were found at 
different projected distances, so that we have in our sample both 
NGC 3073 which is likely strongly affected by its large companion 
galaxy, as well as ID 0650 for which no neighbor is found within a 
0.5° projected radius. Their ellipticities are in the 0.10-0.27 range. 
No kinematic analysis has thus far been available for 2 out of 3 of 
the targets (ID 0650 and ID 0918). 

In total, the chosen galaxies have surface brightnesses at 
the edge of the SAURON field that are not lower than /iy ~ 
24.0mag/arcsec 2 . They are nucleated and non-nucleated objects, 
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Table 1. (1) Galaxy names as in the Virgo Cluster Catalog (VCC) of |Binggeli e taT| {T985), New General Catalog (NGC), or as a GALEX identification number 
(ID), (2,3) J2000 coordinates (from NASA Extragalactic Database (NED), http://ned.ipac.caltech.edu/), (4) for Virgo galaxies: projected distances from the 
cluster's center jLis ker et al. 2 007), (5) surfac e b rightness fluctuatio n (SBF) or globular cluster luminosity function (GCLF) distances for Virgo galaxies where 
avail able jMei et al.|2007| , |Jerjen et al.|2004f or | Jordan et al.|2007f/ ) otherwise an average Virgo distance assumed 4 , for NGC 3073 a SBF distance from |Tonry| 



|et al. 



with 



2001 1, the distances to ID 0650 and ID 0918 were calculated using their radial velocities, all distances are corrected to HO = 73 km/s/Mpc for consistency 
Mei et al.|{2007) , (6) morpholo gical types jLisker et al.|2007| for the Virgo objects, Michielsen et al. 2008 / NED database for the field galaxies), (7) 



ellipticities at 1 effective radius from Lisker et a lj^2007 



R 1 /" growth curves, as detailed in 
on-source exposure times 



Falcon-Barroso et al. 



, (8) r -band effective radii determined by fitting aperture photometry profiles on SDSS images with 
2011 , (9) r-band apparent Vega magnitudes, (10) position angles of our observations, (11) and total 



object a S Rmh7 distance type e R c m r PA t cxp 

(hh mm ss.ss) (dd mm ss.s) (°) (Mpc) (") (mag) (deg) (h) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) 
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Notes: 1 other name: UGC 08986, * other name: CGCG 020-039. 



half of which were known to host a disk/have a substructure. All but 
two have additional spectroscopic (long-slit) and photometric data 
available. The positions, observing details, and photometric proper- 
ties for each observed object are listed in Table[T| 

2.2 Observations 

The observations were carried out over 8 night^] (3 in January 
2010 and 5 in April 2011) with the SAURON integral field unit 
mounted on the 4.2m William Herschel Telescope at the Roque 
de los Muchachos Observatory, Spain. We were able to cover the 
galaxies out to about one effective radius (R e ). The footprints of all 
the pointings overlaid on SDSS images are shown in Figure|2] 

SAURON has a wavelength range of 4760-5300 A. It includes 
a set of potential emission lines (Hj8,[OIII],[NI]) and stellar absorp- 
tion lines (Hj8, Fe5015, and Mgb) which can be used to investi- 
gate stellar populations. We used SAURON's low resolution mode, 
which gives a field of view of 33x41 arcsec 2 , spectral resolution 
(FWHM) of 3.9 A and spatial sampling of 0.94x0.94 arcsec 2 per 
each of the 1431 object lenslets. 

Each galaxy was typically exposed for 10x1800 s. The ex- 
ceptions are 2 field galaxies where 0.5/1. Oh was added due to 
their slightly lower surface brightnesses, and 2 Virgo targets where 
0.5/1. Oh was lost due to problems with rising cryostat tempera- 
ture during our first run. Small offsets were introduced among the 
pointings in order to improve sampling as well as deal with the bad 
columns in the CCD. A dedicated set of 146 sky lenslets, pointing 
1.9 arcmin from the main field, was used to perform accurate sky 
subtraction. As our objects have on average surface brightness a 
few magnitudes below the sky level at the edge of the field, we have 
imaged a number of blank fields ( 1 per night on average) to track 



out of a total of 12 that had been awarded to this project - we lost 3 nights 
in Apr 2010 and 1 in Apr 201 1 due to adverse weather conditions. 



sky variability and optimize sky subtraction. To allow for good cal- 
ibration of line-strength indices we observed a number of stars dur- 
ing each run, covering a broad range of spectral types. Specifically, 
we included stars from the MILES library I Sanchez-Blazqu ez et al.| 
2006) in order to calibrate our measurements. We also observed 
spectrophotometric standard stars each night in order to calibrate 
the response function of the SAURON system. 

2.3 Reduction 

We followed the procedures described in Bacon et al. (200T) for 
the extraction and calibration of the data using the specifically de- 
signed XSAURON software developed at the Centre de Recherche 
Astrophysique de Lyon (CRAL). 

First, the raw frames were preprocessed in a standard way 
which typically includes overscan, bias and dark subtraction. Since 
the dark current values for the CCD used (EEV13) are less than 
1.0 e _ /pixel/h, its effect on the exposures is negligible. 

The next step was to create an extraction mask, i.e. a model 
that relates each pixel on the CCD to its corresponding lenslet and 
wavelength. The resulting table, specifying lenslets' geometrical 
positions on the CCD and their profiles, was used to extract spectra 
from all preprocessed frames and create raw datacubes (a, S, X). 
The cubes were then wavelength calibrated using neon (arc) lamp 
exposures, which contain 1 1 useful emission lines in the SAURON 
wavelength range. These neon exposures, taken before and after ev- 
ery target frame, were used to estimate the possible small offsets 
between the science frame and the extraction mask. This was done 
using a cross-correlation function between its two arc exposures 
and the arc exposure used for the mask creation. 

Flatfielding was a two-step procedure. The regular flatfield 
frames were created one for each observing night using twilight ex- 
posures (spatial flatfield) and continuum (tungsten lamp) exposures 
(spectral flatfield). Since our galaxies fall a few orders of magnitude 
below the sky level at the edge of the field, an additional flatfielding 
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Figure 2. SDSS color images (80x80", or approximately 6.4x6.4 kpc) with the SAURON FoV overplotted. North is up and East left. 



step was included. We took several 900s "blank field" exposures by 
pointing the telescope at relatively empty points in the sky. We then 
created a median of the frames from different days to improve the 
signal-to-noise ratio. The resulting datacube was smoothed at each 
wavelength (i.e. each monochromatic frame) and then divided by 
this wavelength's median to obtain a normalized cube. Because we 
could not exclude night-to-night variability, we decided to check 
how combining certain subsets of sky frames into a superflat and/or 
scaling the gradient in the combined product would affect our sci- 
ence frames. Such an adjustment was needed for our Run 2 objects 
since we opted to use Run 1 skyframes due to their better quality. 
Once the optimal superflat was obtained, we multiplied our regular 
flat fields by the superflat, and these modified flats were incorpo- 
rated back into our pipeline. Dividing the cubes by the superflat 
removed the remaining flat-field residuals. 

Cosmic rays were removed by filtering out statistically impos- 
sible values using a wavelength direction filter run on spectra di- 
vided by the median of a few spatial neighbors. 

Spectral resolution varies slightly within the field of view, with 
the highest values in the center and lowest in the outskirts of the 
field. We therefore measured spectral resolution of each lens by 
comparing our twilight exposures with a high-resolution solar ref- 
erence spectrum and then homogenized it to a common value which 
was less than or equal to the value obtained for 95% of the lenses. 
This was done through applying a gaussian corresponding to the 
difference of the squares of the actual and target resolutions. This 
meant using a value worse than the actual one for some parts of the 
field but it is a necessary step, performed to avoid irregular broad- 
ening of lines when one needs to combine images of an object lo- 
cated at different positions on the chip. We adopted the value of 
FWHM=3.89 A for both runs. 

We performed sky subtraction by computing a median value of 
the 146 dedicated sky lenslets and then subtracting that value from 
the object spectra. The superflatfielding procedure described above 
helped to optimize this step since the sky lenslets spectra (located 
at the edge of the CCD) were the ones most affected by any initially 
present flatfield irregularities. 

Our spectrophotometric standards were used to perform flux 
calibration. A flux correction curve was created to correct for a 
different instrumental response for different wavelengths. This was 
done by comparing the observed spectrum of a flux star with that of 
a reference spectrophotometric spectrum and then applying the re- 



sulting curve to all science frames (see also Kuntschne r et al.|2006| 
for more details). 

We checked the accuracy of our flatfielding process by let- 
ting our twilight frames go through the same reduction steps as 
the galaxy frames (with the exception of sky subtraction since the 
twilights themselves are used for that purpose) and then building 
histograms of the values across the individual frames, which were 
expected to be essentially flat. Indeed, the variations for all tested 
files are always within 1%, and within 0.5% on average. 

To merge the fully calibrated cubes we truncated them to a 
common wavelength range. Then, the cubes' spatial coordinates 
were recentered and the fluxes of corresponding lenslets were co- 
added using a median filter. The resulting cubes were in the process 
resampled to a spatial scale of 0.8x0.8 arcsec 2 . 



3 STELLAR KINEMATICS 

3.1 Stellar kinematic measurements 

In order to accurately measure mean stellar velocity V and veloc- 
ity dispersion cr, a certain minimum signal-to-noise (S/N) level is 
required. To achieve this, we spatially binned our data using the 
Voronoi two-dimensional binning algorithm developed byjCappeF] 
|lari & Copin (2003). In this method one starts from the central bin 
(with the highest S/N) and accretes closest neighboring lenslets un- 
til the required S/N has been reached, then starting a new bin. We 
imposed a minimum S/N of 30 for the kinematic measurements 
presented here. Each time, before binning data data, we also des- 
elected those spatial pixels (spaxels) for which the SfN spaxe i was 
lower than the set minimum (S/N spaxe [ is, effectively, a combina- 
tion of local surface brightness and the number of frames that fall 
on a given spaxel). This extra step was needed since for objects as 
faint as ours, high S/N is critical to recovering accurate kinematic 
parameters values. Simply adding up the signal of individual spax- 
els to form larger bins will lead to lower overall S/N if some of the 
input spaxels have low S/N, and, in our case, showed as a spread 
of (7 values that exceeded the nominal measurement errors. Thus, 
we wanted to know what the minimum SfN spaxe i was that would 
still allow us to recover accurate kinematic measurements. After 
trying out a few values, we adopted SfN spaxe i=7, which roughly 
corresponds to surface brightness of fly ~ 23.5 mag at the edge of 
the fields (slightly depending on the object), as a minimum accept- 
able value. By filtering out low-S/N spaxels we made sure that our 
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Figure 3. Example pPXF fits to the galaxy spectra collapsed within 3 arcsec to illustrate the high quality of the data and the accuracy of the fits. Galaxy spectra 
are shown in black and the fitted models in red. The dark blue lines show the regions used for fitting. The ordinate values are plotted on normalized scales. The 
S/N of these spectra exceeds 100 per pixel. For NGC 3073, emission-line regions were excluded from the fits. In the small panels below each galaxy spectrum 
we plot the residuals between the data and the model, shown as a percentage of the original value (with the y-axes ranges of ±2 %). 
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Figure 5. Stellar velocity dispersion recovery tests. Left panel: the differ- 
ence between the input and recovered o~ is plotted against the input value. 
Right panel: the same difference is now plotted against the signal-to-noise of 
the input spectrum. In both plots the different shades of blue correspond to 
1, 2, and 3 o~ levels of certainty and the dotted vertical lines show the min- 
imum signal-to-noise level, S/N=30, of our bins. The tests in both panels 
only consider the lowest o~ and S/N regimes of our data. As shown in Fig- 
ures|5]and|4] the quality of our data and the comparison with the literature 
are generally much better. 

final binned spectra were not contaminated by low-quality measure- 
ments. 

In the present study, stellar absorption-line kinematics were 
derived for each galaxy by directly fitting the spectra in the pixel 
space using the penalized pixel-fitting method (pPXF) of |CappeL| 
|lari & Emsellem (2004). The method fits a stellar template spec- 
trum convolved with a LOSVD to the observed galaxy spectrum in 
pixel space (logarithmically binned in wavelength). The algorithm 
finds the best fit to the galaxy spectrum by employing nonlinear 
least-squares optimization, and returns the mean velocity V, veloc- 
ity dispersion a, and, if possible, the higher order Gauss-Hermite 
moments and A4, which describe the asymmetric and symmetric 
departures of the LOSVD from a simple gaussian, respectively ( |van| 
|der Marel & Franx|1993| >. 

The procedure of creating the optimal combined stellar tem- 
plate spectrum was repeated for each bin using stellar spectra from 
the updated MILES library (Falcon-Barroso et al. 201 1 1 containing 
spectra of 985 stars at 2.50 A FWHM. From these we selected a 
subsample of ~ 150 stars in such a way that we still had a uniform 
parameter coverage. 

Error estimates were obtained through Monte-Carlo simula- 
tions. For each binned spectrum we ran a loop of 100 realizations, 
each time creating a simulated input spectrum by taking the origi- 
nal spectrum and perturbing it using a spectrum of random values 
drawn from a range determined by the difference between the orig- 
inal data and the best-fit model. Once we had the 100 realizations, 
we built histograms and took the mode with the larger of the error 
bars (since they were asymmetric) as the final computed error value 
for each bin. 

Example fits to the spectra of central regions 1.5 arcsec) of 
all objects are shown in Figure|3] The galaxy spectra are shown to- 
gether with their best-fit models and the data-model residual plots to 
illustrate their excellent agreement and the high quality of the data. 

3.2 Stellar velocity dispersion & comparison with literature 

When the stellar velocity dispersion becomes smaller than the in- 
strumental resolution, a <100 km/s in the case of SAURON, a high 



S/N is required to accurately recover a values. This was extensively 
discussed in, e.g. Toloba et al. « 20 11) who carried out Monte-Carlo 
simulations to determine the minimum S/N needed to recover reli- 
able V and o". Mimicking the instrumental setup of their data, they 
created a set of model galaxies for a range of ages and metallicities 
and found that they were able to recover values with a <10% error 
for all their input a values if the S/N was higher than or equal to 20. 

We repeated the exercise of Toloba et al. (2011 1, and created 
a set of simulations to characterize the uncertainties and biases of 
our velocity dispersion measurements in the low-c and low-S/N 
regimes. A MILES SSP model ( |Vazdekis et al.|2010) of 5 Gyr and 
Solar metallicity (i.e. typical of the average population found in our 
sample) was used as input galaxy. The spectrum was transformed 
to mimic the SAURON instrumental setup of our observations, and 
later convolved to different velocity dispersions and signal-to-noise 
ratios. The results of this experiment are shown in Fig. [5] In the left 
panel we show the potential biases in the recovery of the velocity 
dispersion for the lowest S/N considered in this work. As expected, 
the average difference between the input and output values tends 
to zero as o"i nput increases. For o"j nput below our pixel scale (~60 
km/s) one can start to appreciate systematic behaviors (i.e. over- 
and under-estimation of the input values). This effect is, however, 
well below the uncertainties at the lowest velocity dispersions we 
measure in our maps (-~30 km/s). The right panel shows typical un- 
certainties for a large range in S/N ratios in the velocity dispersion 
range where the line-of-sight velocity distributions are undersam- 
pled. As in the previous plot, at the minimum signal-to-noise level 
considered for the Voronoi binning any systematic trend is within 
our measured uncertainties. These tests suggest that our measure- 
ments should not be significantly affected by sampling effects. 



To directly establish the level of confidence of our results 
we compared our a measurements with literature values. Figure [4] 
shows the comparison between our and previously published pro- 



files of Simien & Prugniel ( 2002 1, Geha et al. (|2003|l,|van Zee et al. 


2004 


i,|Beasley et al.|(|2009f, Chil ngarian (2009), and Toloba et al. 


2011 


1, typically based on high spectral resolution measurements. 



We plot both the individual bin values, as well as values calculated 
by averaging the bin values over ellipses (using the kinemetry soft- 
ware of Krajnovic et al. 2006). The errors shown are representa- 
tive of average bin errors at a given radius. For most of our objects 
the agreements is good and the profiles agree. For VCC 1036 and 
VCC 2048 the values generally agree at larger radii but we do not 
see a drop in o" reported by some authors. We do, however, agree 
to within the errors with the published profiles of VCC 1261, for 
which a a drop has also been reported. Certain disagreement be- 
tween our and published data for VCC 0929 is very likely caused 
by the lower quality of the published data, evidenced by large scat- 
ter and error bars of those measurements. The overall conclusion 
from the comparison is that we are able to measure a values which 
agree to within the errors with most of the profiles available in the 
literature. Thus, the claim that an accurate determination of <T as 
low as ~half of the instrumental resolution is possible with high- 
quality data is now supported not only from a theoretical but also 
an observational point of view. 



4 LINE STRENGTHS 

4.1 Line strengths measurements 

The measurement of absorption line strengths in combination with 
stellar population models has been used for years to probe the 
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Figure 4. Comparison of the published o~ profiles with the ones obtained in this study. The light-gray circles show all our bins and in red we show profiles 
obtained by averaging the bin values over ellipses, the overplotted errors are average errors for bins at a given elliptical distance. Literature values are shown 
using different colors/symbols as denoted in the VCC1861 panel. The following abbreviations were us ed: l)|Pedraz et al.|(2002) : P02, 2)|Simien & P rugniel 
j2002) : S02, 3) |Geha et al.lpOtB) : G03, 4) |van Zee et"aI]j2fJ04) : Z04, 5) |Beasley et"aH{2009) : B09, 6) |Chilingarian|p009| : C09, and 7) |Toloba et al.|(201 1) : 
Tl 1 . For VCC 0929, VCC 1036, and VCC 1261 some of the published profiles extend beyond what is presented in the figure. 
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Figure 6. The comparison between the line-strength values obtained for the MILES stars we observed and the values for the same stars from the MILES stellar 
library of Sanchez-Blazquez et al. 1 2006). The thick dashed red lines show the mean value of the difference between the library and our values, and the thin 
dotted-dashed red lines indicate the standard deviation. Zero is shown with the dotted black line. All values are in LIS-5.0 A system of |Vazdekis et al.|{2010} . 
See the main text for more details. 



luminosity-weighted age, metallicity and abundance ratios of cer- 
tain elements in integrated stellar populations. 

For many years now the standard in the field has been the set 
of 25 indices in the Lick/IDS system fWorthey et al.|1994fr . Unfor- 
tunately, the system has been defined on a low-resolution library 
(8-11 A), which with the advent of high-resolution spectroscopic 
data meant that in order to use the indices one had to degrade the 
resolution of their data to match that of the Lick/IDS system, thus 
effectively forsaking the data quality. Therefore, a need emerged to 
define a new system that would avoid this data degradation, as well 
as other problems inherent to the Lick library, such as changing 
resolution as a function of wavelength, uncertain calibration of the 
continuum shape, and low effective S/N due to significant flat-field 
noise (see e.g. Worthey et al. 1994 or Trager et al. 1998). 



The new Line Index System (Vazd ekis et al.||2?)IO) meets 
these requirements. It provides constant resolution over the whole 
wavelength range, and it avoids the above mentioned intrinsic 
uncertainties associated with the Lick/IDS system. The system is 
defined at three different resolutions, 5.0 A (suitable for globular 
clusters' index comparisons), 8.4 A (appropriate for studies of 
dwarf and medium-mass galaxies, roughly matches the Lick/IDS 
system resolution), and 14.0 A (designed for studies of massive 
galaxies). The MILES models are based on the MILES stellar 
library presented in Sanchez-Blazquez et al. (2006) that contains 
985 stars covering a much wider parameter space (T e ff> l°g 
(g), [Fe/H]) - particularly at lower metallicities - than any other 
library available before. Their sp ectral range is 3540-7610 A at 
2.50 A JFalcon-Barroso et al.|201 1 1. 



What indices can be measured is mainly determined by the 
available wavelength range of one's instrument and a given object's 
redshift. When observing Virgo objects with SAURON we have 3 
traditional Lick indices available: Hj3, Mgb and Fe5015, of which 
Hj3 is the age-sensitive index and the other two are metallicity in- 
dicators. In this paper we present line-strength maps and profiles in 
the LIS-5.0 A system. 

The errors on the indices were calculated through Monte-Carlo 
simulations. For each realization, we perturbed the spectra to incor- 
porate the following error sources: V measurement, which affects 
the position of the lines; o measurement, which affects the lines' 
depth; and systematic errors (in the same way as described earlier 
for kinematic measurements, using as the error spectrum the differ- 
ence between the data and the best-fit model spectra). 



4.2 Comparison with published values 

Central integrated line-strength values have previously been pub- 
lished for a number of dEs by|Geha et al.] p003 1, Mic hielsen et al.| 
p008") , and |Paudel eTaT1 ( |2010) . While we do have a few galax- 
ies in common with their samples, a comparison between their and 
our line-strength values is not straightforward. The published values 
are in the old Lick/IDS system and the conversion from Lick/IDS 
to LIS is rather uncertain and may introduce biases. Moreover, to 
properly convert we would need to know the exact instrumental se- 
tups used to obtain the published values. 

We thus opted for a different method to test how reliable our 
line-strength measurements are. We measured the indices on the 
MILES stars we observed every night and then compared their val- 
ues with those of the same stars in the MILES stellar library of 
Sanchez-Blazquez et al. ( 2006 ), measured at our instrumental reso- 
lution using the MILES website (miles.iac.es). By choosing to test 
our measurements in this way, we made sure that we had control 
over all factors to which line-strength measurements are sensitive: 
the overall data quality, reduction procedures, spectral resolution, 
and, most importantly, the system in which line strengths were mea- 
sured, without the need to resort to conversion tables. 

The results of the comparison are in Figure [6] where we show 
the absolute differences between the observed and the library stars 
values. The agreement for Hj3 and Mgb measurements is very 
good, with the average offsets being only 0.005±0.110A and - 
0.03±0.147 A, respectively. The offset for Fe5015, though a littler 
higher (0.189±0.370 A) is still, like in the case the other two in- 
dices, within the error bars of our galaxy observations (see the ap- 
pendices). The means and standard deviations of the relative errors 
((indexoBS — m dexMiLEs)/ m dexiy[iLES> n °t shown in the figure) are 
0.003±0.056, 0.063±0.108, and -0.020±0.054 for H/3, Fe5015, 
and Mgfo, respectively. We thus conclude that our measurements 
are robust. 



5 RESULTS 

With integral-field data we are able to build not just kinematic (V 
and c) and line-strength profiles, but maps that provide information 
for every spatially resolved element. This is not the first study de- 
voted to dEs in which integral-field data has been used (see Prugniel 
|et al.||2005| |Chilingarian |2009> , however, our maps cover a much 
larger area than those few published previously so we can charac- 
terize the behavior of kinematic and stellar population parameters 
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Figure 7. Stellar velocity and velocity dispersion maps of our full sample. The maps are overplotted with the isophotes of the reconstructed images and 
presented on a common scale to facilitate the comparison among galaxies. The same maps (including errors) are presented in Appendix [A] with optimized 
scales to emphasize the individual characteristics. 



out to much (up to four times) larger radii. The available long-slit 
data do sometimes go a bit farther out ( Chilingarian 2009 Toloba| 
|et al.pOTT) but the data never reach beyond 25-30". The spatial 
information obtained from our data is, thus, much more extensive. 
Also, when analyzing radial properties, we can average our data 
points over annuli, which significantly increases S/N and hence re- 
duces error bars. 

Appendix |A| contains a full set of kinematic and line-strength 
maps, including error maps, for each galaxy. Additionally, in Ap- 
pendix [B] we plot and tabulate one-dimensional profiles of all mea- 
sured quantities and their errors in order to enable an easier compar- 
ison of our and literature results. Here we show all maps together 
to emphasize the overall variety and to facilitate the comparison of 
properties among the galaxies. The rest of this section provides a 
detailed description of these properties. 



5.1 Stellar kinematic maps 

Figure [7] shows our V and a maps. The scale is the same for 
all galaxies to emphasize their diversity. The maps show vary- 
ing degrees of rotation: from V max ~40 km/s for VCC0523 to 



virtually no rotation for VCC1261. These rotation values seem 
to be uncorrelated with flattening since we see all possible com- 
binations of both: flattened rotators, flattened non-rotators, fairly 
round rotators and, similarly shaped non-rotators. In addition to 
that, in two galaxies (ID 0650 and ID 0918) we observe large-scale 
kinematically-decoupled components, counter-rotating with regard 
to the main bodies. The flattened slow rotator VCC 1087 and non- 
rotator VCC 1261 are known (through globular clusters analysis 
of |Beasley et al.|[2009) to display strong rotation at much larger 
radii (4-7 R e ). The rotation in the (almost) round VCC 0308 and 
VCC 0929 could be due to them being seen not perfectly face-on, 
or it could be ascribed to their specific orbital configuration if they 
were triaxial. We see kinematic twists in two systems: in VCC 0523, 
known to h ost a bar (TLisker et al.|2 007 ) and showing weak signs of 
spiral arms (Lis ker et al.|2009 1, it is very clear that the central pho- 
tometric isophotes are twisted with regard to the outer ones, and 
it is the outer ones that agree with the kinematic major axis of the 
galaxy; in ID 0918 the twist is less pronounced, but has similar char- 
acteristics. 

The velocity dispersion maps are rather flat, in the o~ ~40-60 
km/s range, with shallow gradients for a few of the galaxies: a slight 
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central increase (VCC 1431) or decrease (VCC0929, VCC 1261). 
Only for 1 field galaxy is a significantly higher (ID 0918, ~80 
km/s). The interpretation of flat dispersion profiles requires employ- 
ing dynamical models and will be presented in our Paper II. 

What is most significant in these stellar kinematic maps is that 
no two galaxies in our sample are really alike. The implication of 
this to their formation histories is discussed in Section|6] 



5.2 Line-strength maps 

The goal of the stellar population analysis is to gain insight into the 
star-formation and assembly history. For example, two-dimensional 
data makes it possible to connect stellar populations with kinematic 
structure, as shown by Kuntschner et al. (2006 ) where various dis- 
crepancies between isophotal and isoindex contours were found for 
their sample of giant elliptical (E) and lenticular (SO) galaxies. 

We do not seem to find such strong signatures in our data 
(Figure [8}, i.e. our data are consistent with no deviation between 
isophotes and isoindex contours with the exception of one galaxy 
(VCC 1036, Fe5015 map). This would suggest that the fast-rotating 
components, i.e. disks, are not different from the main bodies in 
terms of populations. However, the gradients in the indices do dif- 
fer among our galaxies: for Fe5015 and Mgb (i.e. the metallicity 
tracers), while negative (i.e. more metal-rich in the center) in all 
galaxies but NGC 3073 , they do vary from almost flat to having 
a strong peak (or in the case of NGC 3073: a dip) in the center; 
Hj3 (age tracer) gradients are similar in that respect, some show- 
ing a highly-concentrated peak (like in the blue-core VCC 0308), 
with others being essentially flat. The general trend we get from our 
maps is that dEs are on average older and more metal-poor in their 
outer parts, in agreement with Michiel sen et al.|j2008} . We further 
investigate these points in Paper III, where a detailed stellar pop- 
ulation analysis, including population gradients and star formation 
histories (using full spectrum fitting techniques) will be presented, 
and in Paper IV, where we use aperture measurements to analyze 
global scaling relations and compare dEs to other galaxy types. 



5.3 Notes on individual galaxies 

Here we provide a brief description of the properties of our maps 
for each galaxy. The following abbreviations will be used through- 
out the section: 1) Pedraz et al. (2002): P02, 2) Simien & Prugniel 
pool} : S02, 3) |Geha et al.| l |2003) : G03, 4) |vanZee et al.H2004) : 
Z 04, 5)|Lisker et al.|p006l >: L06, 6) |Chilingaria"n] ( (2009l ): C09, and 
7) |Toloba et al.H2011) : Tl 1. 

VCC 0308 is designated dE(bc,di) in the notation of L06, and 
it is the only blue-core (be) galaxy in our sample. The Hj8 map 
shows a strong peak in the center, suggestive of the presence of a 
younger population, consistent with the '(be)' classification. It is a 
round galaxy in the outskirts of the cluster, showing a mild level of 
rotation, signs of which have also been seen by Tl 1. 

VCC 0523 is a mildly flattened, nucleated dwarf in the inner 
parts of the cluster; L06 found a disk and a bar in this galaxy based 
on photometric decomposition, so the galaxy is denoted dE(di). We 
confirm this result by showing how the photometric and the kine- 
matic major axes of the galaxy are misaligned due to the bar. The 
galaxy is rotationally supported, in agreement with Tl 1 results from 
their mean velocity and velocity dispersion curves. 

VCC 0929 is a fairly round galaxy with a high level of rota- 
tion. Previous results also indicate rotation, though to a lesser de- 
gree (P02) or showing large uncertainties (S02). Line-strength maps 



are essentially flat, showing only a mild increase towards the center 
in the Fe5015 and Mgb maps. 

VCC 1036, of the dE(di) type, is a very flat object that lies in 
the central parts of the cluster. It is a well-studied galaxy, previously 
analyzed by S02, G03, Z04, and C09. In agreement with these stud- 
ies we find that the galaxy is fast-rotating. We do not confirm the 
presence of a central a drop reported by C09. The HjS map is flat, 
and Fe5015 map shows a central increase/peak, which is yet more 
pronounced in the Mgb map. 

VCC 1087, centrally located, is another well-known flattened 
nucleated dwarf, this time, however, with no substructures (L06) 
and little apparent rotation as reported previously by S02, G03, C09 
and Til. Our maps do show some level of rotation at the edge of 
the field, it is, however, almost consistent with zero. Interestingly, 
Beasley et al. (2006) results from the analysis of globular clusters 
located at up to 4-7 effective radii reveal a significant level of rota- 
tion at those larger distances. It is, thus, possible that our maps look 
flat due to insufficient radial coverage. 

VCC 1261 is the only galaxy in our sample, which, despite 
significant flattening, truly exhibits no signs of rotation out to 1 
effective radius (previously noted by G03, Z04, C09 and Til). 
C09 velocity curve and sigma profile suggest a possible presence 
of a small kinematically decoupled component. Our data, as 
well as those of G03 and Til, do show some signs of a central 
<T drop, though it is more spatially extended and not as deep. Our 
line-strength maps show a central peak in the H/3 and Mgb maps, 
suggestive of a rather strong positive age and negative Z gradient. 
The galaxy's GCs were studied by Beasley et al. (2009) with the 
similar findings as for VCC 1087 above. 

VCC 1431 is a very round, non-rotating galaxy in the center 
of the cluster. The galaxy has previously been analyzed by Til. 
We note a small increase in cr towards the galaxy's center, which is 
however, consistent with zero to within the errors. The line-strength 
maps show a clear central peak in the Fe5015 map, while the Hj8 
and Mgb maps are flat. 

VCC 1861 is an almost perfectly spherical non-rotator (with 
an average V < 10 km/s), but it must be remembered that, like in the 
case of VCC 1431 above, we may be viewing the galaxy face-on, 
in which case we would be unable to detect any (possible) rotation. 
No substructures are present as reported by L06. Our line-strength 
maps show a small negative Fe5015 and Mgb gradient. 

VCC 2048 is a flattened disky rotator situated in the clus- 
ter's outskirts. Previous kinematic analyses are found in S02 and 
C09. The latter finds a pronounced a drop in the center. The drop 
would be undetectable for us at the claimed minimum cr level 
(20 km/s). Our line-strength maps seem fairly flat, with the excep- 
tion of the MgZ? map which shows a slight increase towards the 
center. VCC 2048 is the galaxy considered by Kormendy & Bender 
(201 1) to be a transition object and a "missing link" between the SO 
and dE (Sph in their notation) classes. 

NGC 3073 is a field galaxy in the vicinity of the larger 
NGC 3079 (with their projected separation of ~10 arcmin). The 
galaxy has previously been observed by Til and Emsellem et al. 
( |201 1| >. The classification of this galaxy varies from source to 
source: SAB0~ in |de Vaucouleurs et al.|{l99~i) , SO in |Nilson|{l973"t , 
and Til classify it as dE(di). The galaxy's integrated intensity 
isophotes are distorted by the presence of gas which is later seen in 
the spectra in the form of strong emission lines . |Serra et aL] ( |201 1) 
show the Hi pattern for NGC 3073 where a tail is present pointing 
away from its large neighbor. Small level of rotation was previously 
detected by S02 and Til, and is indeed confirmed here. 
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Figure 8. Line-strength maps of our full sample. The maps are overplotted with the isophotes of the reconstructed images and presented on a common scale to 
facilitate the comparison among galaxies. See AppendixlAlfor maps with individually optimized scales and errors. 



ID 0650 is a truly isolated field galaxy, with no galactic neigh- 
bors within ~30 arcmin projected radius. This is the first kine- 
matic data for this object available in the literature. We find a large- 
scale counter-rotating component with an approximate diameter of 
12arcsec (or ~1.5 kpc). The Fe5015 and Mgb maps show a slight 
central peak. 

ID 0918 is a field galaxy with the closest galactic companion 
at a projected distance of ~20 arcmin, also with no kinematic data 
available so far. The integrated intensity isophotes show a slight 



misalignment of the inner and outer parts, the latter agree with 
the galaxy's rotation pattern. We report a discovery of a strong 
counter-rotating component, whose approximate size is 14arcsec 
(or ~ 1.1 kpc) along the major axis. The line-strength maps show a 
strong central peak in the Fe5015 and Mgb maps. 
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Figure 9. Integrated (within 1 effective radius) specific angular momentum 
Xr versus projected (open black circles) and deprojected (filled blue circles) 
Virgocentric distance. The errors on A/j e are assumed to be 0.05 after |En>] 
sellem et al. (2011 ) (see their Section 4.5). The dotted lines show the fits to 
the data for both types of distances. While a weak trend is seen when the 
projected distances are used, this is not true for the intrinsic deprojected dis- 
tance. No relation between the rotational support and location in the cluster 
is thus seen in our data. 



6 DISCUSSION AND CONCLUSIONS 

We have presented stellar kinematic and absorption line-strength 
maps for a sample of 12 Virgo Cluster and field dwarf elliptical 
galaxies. This is to date the largest sample of dEs for which integral- 
field data is available, and offers the largest spatial coverage so far. 
We have shown how having high S/N data and performing care- 
ful data reduction enables one to recover velocity dispersion values 
well below instrumental resolution. We do not discuss here all the 
relations presented in the Introduction since this will be done in 
more detailed subsequent papers that will each focus on different 
aspects of dEs. Instead, we outline some of the most noteworthy as- 
pects of the data and focus on the meaning of the objects' structural 
variety. 

Even having a relatively small sample we still find all possible 
combinations of flattening and rotation in our kinematic maps, to- 
gether with twists and decoupled components. The maps not only 
show an overall remarkable variety of properties but also draw our 
attention to interesting individual cases. Our flattened slow/non- 
rotators (VCC 1087 and VCC 1261) deserve particular attention 
since the analysis of their globular clusters systems (Beasley et al. 
|2006| and |Beasley et al.|2009| l shows rapid rotation at larger radii 
and, in one case, a slight misalignment of axes between the GC and 
steflar components. One possible explanation could be that the ini- 
tially present angular momentum has migrated outwards, e.g. due 
to a disk-halo coupling, still, the actual processes are difficult to 
name. 

The kinematicaily-decoupled components (KDCs) found in 
ID 0650 and ID 09 18 are similar in physical size to those found in 
giant elliptical galaxies of Krajn ovTc et al.|p01 1} . KDCs in dwarfs 
are more difficuft to explain given that we generally believe they 



have not formed through mergers (De Rijcke et al. 2004). For clus- 
ter members there exist studies that suggest the possibility of cre- 
ating a KDC via harassment ( Gonzalez-Ga rcfa et al.| 2005 I. Still, 
for isolated galaxies like ID 0650 and ID 09 18 the puzzle remains 
unsolved, since ram-pressure stripping, though believed to be able 
to act even on objects outside of dense environments (especially if 
concepts like warm-hot intergalactic medium are invoked, e.g. [Dave] 
|et al.|200"T| which would allow for more efficient stripping), is nev- 
ertheless unable to sufficiently alter the stellar composition. It is yet 
unclear whether the peaks in the Fe5015 and Mgfe maps are related 
to the presence of the KDCs, i.e. whether steilar-populations-wise 
the KDCs could be distinct from the main bodies of the galaxies 
(e.g. external gas accretion can be responsible for creating counter- 
rotating components). Ail these findings invite further investigation, 
in particular of the orbitaf make-up of these systems. 

The line-strength maps also show a significant degree of di- 
versity. The gradients visible in the maps vary from flat to strongly 
peaked in the centers, showing a spread in stellar ages and metallic- 
ities, both among and within the galaxies. All but one systems are 
essentially gas-free, indicative of no recent or ongoing star forma- 
tion (the exception is one of the field galaxies, NGC 3073). Regard- 
less of the variety, and as has already b een found by |Michielsen| 
|et al.| (|2008), the comparison of our and Kunts chner et al.| ( |2010} 
integrated line-strength values shows that the ages of dEs are on av- 
erage younger, and the metallicities lower than those of giant early- 
type galaxies (ETGs, Rys & Falcon-Barroso 2012). What is worth 
noting is that, unlike for a large number of giant early-type galax- 
ies of Kuntschne r et al.| (2006 1, no deviations are found between 
isophotes and isoindex contours of our dwarfs. This could suggest 
that any kinematic structures have not been built up by star forma- 
tion activity restricted to a certain region of a given galaxy. We will 
explore these questions in more detail in a later paper of this series. 

To find out whether the presence of photometric and kinematic 
(sub)structures is correlated in any way, we have made one-to-one 
comparisons of our kinematics with the results of Janz et al. (2012b, 
in prep.). We find that among multiple component objects we have 
both rotators and non-rotators, the same goes for single-component 
galaxies; also, two of our slow-rotators contain lenses. Additionally, 
four out of five rotators are disky galaxies, based on the Lisker et al. 
< |2007) > nomenclature. It thus seems that the presence of photometric 
and kinematic (sub)structures is not correlated in any obvious way, 
which yet again strengthens the variety argument. 

We note here that this variety combined with the variety found 
by |Krajnovic et al.| pOTT) in the ATLAS3D ETGs should not be 
viewed as a sign of dEs and giant ellipticals (Es) being continu- 
ous in their properties. Namely, 70% of the galaxies included in 
the ATLAS3D sample are SO galaxies and they have been shown 
to be fundamentally different from the other ETG type, Es, in both 
their kinematic (ATLAS3D and SAURON surveys) and photomet- 
ric (Kormendy & Bender 2012) properties, with SOs (as defined 
photometrically) and fast rotators (as defined kinema tic ally j^] hav- 
ing properties more similar to those of late-type galaxies. If indeed 
SOs are faded spirals, as suggested most recently by Kormendy & 
Bender ( 2012 ), then the fact that our dEs show a similar level of va- 
riety as the ATLAS3D fast rotators only supports the claim of them 
being transformed late-type galaxies. 

Our results are, thus, best understood in the context of forma- 
tion scenarios that, first, assume an infall of younger, late-type ob- 
jects as the source of dE progenitor galaxies, and second, allow for 



The terms are not interchangeable, as discussed in Emsellem et al. 201 1 . 
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the random shaping of galaxy properties, based on the variety initial 
infall parameters (disk inclination and orbit type) and the number 
of encounters with large galaxies. The combined influence of ram- 
pressure stripping and harassment fulfills these requirements, still, 
the exact impact of the two effects is not yet known (see, e.g. the 
simulations of |Smith et al.|20l"2| which suggest that the removal of 
the gas potential could be the key source of disk heating). The tell- 
tale signs of the environmental influence can be found in changing 
galaxy properties as a function of clustrocentric distance, as dis- 
cussed in the Introduction. 

What we would like to stress, however, is that the analysis of 
any distance trend needs to be performed with caution, in partic- 
ular, the depth of a cluster needs to be taken into account since 
using projected (instead of physical) distances might yield biased 
results, especially when dealing with incomplete samples. For ex- 
ample, we do not find a confirmation for the |Toloba et al.| ( |2009} 
claim that rotationally-supported galaxies are preferentially located 
in the cluster outskirts. To determine the level of rotational sup- 
port we u sed the specific angular momentum Xr c of Emsellem et al. 
(2007^ where Re signifies values integrated within 1 effective ra- 
dius, and then plotted it as a function of deprojected clustrocentric 
distance. The relation we find is essentially flat. The same plot using 
projected distances (not shown here) shows a slight upward trend, 
i.e. would be in agreement with Toloba et al. (2009). Still, our sam- 
ple size prohibits any firm conclusions. 

It is also important to note that, while we expect to find all 
these trends relations of properties with distance, we also expect to 
see scatter around them, including counter examples. This is be- 
cause transformation by repeated encounters is, to a certain degree, 
a random process, therefore there will always be galaxies whose 
characteristics have been greatly altered in an individual encounter 
of a specific configuration that, for example, unexpectedly creates a 
non-rotating galaxy in the outskirts of a cluster, a rotating one close 
to the center, a fairly round, rotating object, etc. 

Our findings fit this picture remarkably well. Unlike other 
autho rs (e.g. |Lisker et al.|[2008l |Koleva et al.p009l |Paudel et al.| 
|2011) we believe that the presence of a (large) variety of properties 
does not necessarily imply different formation channels. Since both 
ram-pressure stripping and harassment are stochastic processes, 
a large spread in properties is even expected. Thus, we consider 
the variety of characteristics we find in this work to be our most 
important result. 

The detailed analysis of relations of dEs to other galaxy types, 
their dynamical structure, and their star-formation histories is an 
attempt at providing answers to the question of the relative impor- 
tance of various proposed transformation mechanisms and will be 
presented in the forthcoming papers. It would, in addition, be very 
valuable to study in two dimensions the properties of larger sam- 
ples of irregulars, dwarf spirals, and field dEs, to compare them 
with cluster dEs. 
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APPENDIX A: MAPS 

Here we provide maps of intensity, V, o", and line strengths for the 3 
indices, Hj3, Fe5015, and Mgb, in the LIS-5.0 A system. Each map 
is accompanied by an error map, obtained through Monte-Carlo 
simulations as described in section [3~T] The plots are complemented 
with basic galaxy characteristics, repeated here from Table 1 for the 
reader's convenience. The plotting ranges are adjusted to emphasize 
individual features. 



APPENDIX B: PROFILES 

In order to facilitate the comparison between our and literature data, 
we plot here (Figure [Bl} one-dimensional (ID) profiles of all our 
measured quantities, mimicking long-slit data. Tabulated values to- 
gether with errors are provided in Table [Bl] The profiles were ex- 
tracted along ellipses in steps of 1" using the kinemetry software of 
|Krajnovic et al.| |2006 l. The listed errors were obtained by running 
the software on the error maps presented in Appendix |A| so they 
reflect average bin errors at a given elliptical distance. 
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Figure Al. Maps of intensity, V, o~, and line strengths for the 3 indices: H/3, Fe5015, and Mgfo, accompanied by error maps. See Appendix [Alfor more details. 
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Figure Al. Continued. 
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Figure Al. Continued. 
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Figure Al. Continued. 
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Figure Bl. One-dimensional (ID) profiles extracted from our maps, shown together to illustrate differences among our galaxies and to facilitate the comparison 
between our and literature data (the vast majority of the latter available in ID). The errors, not shown here for clarity, are listed in Table |Bl| In the Hj3 panel 
the NGC3073 profile (with values in the 5.26-3.02 range) is not included since the y-axis range is restricted to emphasize other galaxies' profile shapes. 
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Table Bl. One-dimensional profiles of the stellar velocity V, velocity disperson o~, and Hp, Fe5015, and Mgb lines, and their errors, as a function of elliptical 
radius. The profiles were extracted in steps of 1" using the kinemetry software of Krajnovic et al. 1 2006). 
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1.98 


0.59 


17. 


2.36 


6.5 


52.1 


10.5 


2.00 


0.57 


2.99 


1.28 


1.93 


0.61 


18. 


1.33 


6.1 


53.4 


10.7 


2.07 


0.56 


2.87 


1.27 


1.89 


0.60 


19. 


0.10 


6.2 


56.0 


11.6 


2.05 


0.56 


2.68 


1.28 


1.86 


0.61 


VCC 1431 


0.0 


0.00 


1.0 


45.9 


2.44 


1.91 


0.10 


3.46 


0.22 


2.63 


0.10 


1.0 


4.92 


1.2 


53.2 


3.11 


1.84 


0.14 


3.55 


0.31 


2.80 


0.15 


2.0 


0.75 


1.6 


56.2 


3.79 


1.77 


0.19 


3.64 


0.40 


2.85 


0.20 


3.0 


3.31 


2.2 


55.1 


4.90 


1.72 


0.24 


3.51 


0.52 


2.84 


0.26 


4.0 


2.78 


2.6 


53.1 


6.07 


1.71 


0.28 


3.42 


0.60 


2.76 


0.30 


5.0 


1.23 


3.1 


52.9 


6.81 


1.70 


0.32 


3.52 


0.69 


2.65 


0.33 


6.0 


2.32 


3.9 


52.9 


8.07 


1.65 


0.40 


3.33 


0.84 


2.51 


0.41 


7.0 


2.78 


4.7 


49.7 


8.81 


1.62 


0.46 


3.23 


0.98 


2.49 


0.47 


8.0 


6.16 


5.4 


47.4 


9.32 


1.68 


0.53 


3.20 


1.11 


2.44 


0.55 


9.0 


3.77 


5.7 


48.6 


10.2 


1.80 


0.55 


3.32 


1.15 


2.32 


0.58 


10. 


5.01 


6.3 


47.1 


10.5 


1.69 


0.60 


3.21 


1.25 


2.17 


0.62 


11. 


2.47 


7.0 


45.0 


10.4 


1.69 


0.67 


3.18 


1.37 


2.12 


0.67 


12. 


1.95 


7.6 


43.8 


10.5 


1.68 


0.75 


3.38 


1.54 


2.09 


0.75 
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Table Bl. Continued. 



r 


V 


eV 


a 


eo" 


H/3 


eH/3 


Fe5015 


eFe5()15 


Mgb 


eMgb 


(") 


(km/s) 


(km/s) 


(km/s) 


(km/s) 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


VCC1861 


0.0 


0.00 


1.6 


42.8 


4.02 


2.01 


0.19 


3.83 


0.38 


2.54 


0.20 


1.0 


6.08 


2.4 


38.2 


4.24 


2.00 


0.24 


4.01 


0.52 


2.56 


0.26 


2.0 


1.18 


2.7 


38.8 


4.74 


2.05 


0.28 


4.13 


0.61 


2.69 


0.31 


3.0 


2.84 


3.0 


35.8 


4.95 


2.09 


0.33 


4.12 


0.73 


2.68 


0.35 


4.0 


3.32 


3.2 


36.7 


5.54 


2.11 


0.35 


4.24 


0.76 


2.57 


0.39 


5.0 


4.88 


3.6 


36.5 


5.81 


2.06 


0.40 


4.11 


0.85 


2.54 


0.44 


6.0 


4.23 


4.1 


38.3 


6.49 


2.03 


0.44 


3.80 


0.96 


2.54 


0.48 


7.0 


4.90 


4.3 


37.9 


6.68 


2.09 


0.46 


3.58 


0.99 


2.48 


0.50 


8.0 


6.11 


4.7 


36.9 


7.00 


2.01 


0.50 


3.49 


1.05 


2.42 


0.55 


9.0 


3.96 


4.9 


38.2 


7.41 


2.07 


0.50 


3.41 


1.12 


2.39 


0.57 


10. 


5.10 


5.3 


37.3 


7.61 


2.07 


0.54 


3.27 


1.21 


2.39 


0.60 


11. 


5.44 


5.7 


37.4 


7.70 


2.02 


0.57 


3.16 


1.27 


2.40 


0.62 


12. 


5.19 


5.6 


36.5 


7.50 


2.06 


0.57 


3.18 


1.26 


2.34 


0.63 


13. 


5.81 


5.9 


37.1 


7.96 


2.08 


0.60 


3.13 


1.30 


2.32 


0.65 


VCC2048 


0.0 


0.00 


1.0 


46.5 


2.15 


2.26 


0.10 


3.97 


0.25 


2.02 


0.12 


1.0 


9.55 


1.2 


42.6 


2.85 


2.27 


0.13 


3.50 


0.28 


2.00 


0.14 


2.0 


8.67 


1.5 


41.8 


3.39 


2.36 


0.15 


3.45 


0.33 


1.93 


0.16 


3.0 


5.48 


1.8 


43.2 


3.92 


2.32 


0.18 


3.35 


0.39 


1.88 


0.19 


4.0 


7.91 


2.1 


42.6 


4.46 


2.25 


0.20 


3.29 


0.45 


1.81 


0.22 


5.0 


12.8 


2.3 


45.8 


5.10 


2.25 


0.23 


3.31 


0.49 


1.85 


0.24 


6.0 


15.5 


2.7 


45.6 


5.77 


2.24 


0.25 


3.21 


0.55 


1.80 


0.27 


7.0 


16.0 


2.9 


46.7 


6.13 


2.19 


0.28 


3.15 


0.62 


1.76 


0.31 


8.0 


15.7 


3.3 


48.7 


6.90 


2.15 


0.32 


3.06 


0.71 


1.76 


0.34 


9.0 


13.9 


3.9 


44.2 


7.25 


2.13 


0.35 


2.98 


0.77 


1.76 


0.38 


10. 


13.9 


4.4 


45.6 


8.19 


2.08 


0.41 


3.14 


0.89 


1.71 


0.44 


11. 


16.5 


4.8 


48.4 


9.06 


2.07 


0.45 


3.19 


1.00 


1.67 


0.49 


12. 


13.8 


5.2 


49.8 


9.31 


2.12 


0.50 


3.07 


1.07 


1.73 


0.53 


13. 


17.0 


5.5 


51.4 


10.2 


2.10 


0.53 


2.97 


1.12 


1.75 


0.55 


14. 


21.8 


5.7 


48.1 


9.81 


2.14 


0.55 


2.97 


1.18 


1.69 


0.58 


15. 


20.4 


5.9 


46.5 


9.89 


2.06 


0.58 


2.96 


1.25 


1.68 


0.62 


16. 


20.5 


6.1 


49.4 


10.7 


1.98 


0.62 


2.90 


1.31 


1.72 


0.66 


17. 


20.5 


6.6 


52.5 


12.1 


1.94 


0.67 


2.83 


1.42 


1.72 


0.71 


18. 


17.5 


6.9 


54.2 


13.2 


1.96 


0.73 


3.01 


1.54 


1.81 


0.77 


NGC3073 


0.0 


0.00 


1.3 


49.2 


3.36 


5.26 


0.08 


2.43 


0.20 


0.86 


0.09 


1.0 


0.19 


1.5 


45.5 


3.79 


5.12 


0.09 


2.49 


0.20 


0.99 


0.10 


2.0 


1.33 


1.7 


43.5 


4.18 


4.88 


0.10 


2.45 


0.24 


1.13 


0.12 


3.0 


2.17 


2.1 


46.2 


4.81 


4.66 


0.14 


2.62 


0.30 


1.26 


0.15 


4.0 


3.13 


2.6 


49.0 


5.69 


4.43 


0.18 


2.80 


0.39 


1.36 


0.20 


5.0 


7.82 


2.8 


46.6 


6.26 


4.25 


0.21 


2.85 


0.46 


1.45 


0.24 


6.0 


7.95 


3.3 


45.8 


6.46 


4.09 


0.24 


3.10 


0.53 


1.50 


0.27 


7.0 


5.44 


3.6 


48.1 


7.69 


3.91 


0.28 


3.12 


0.61 


1.59 


0.31 


8.0 


4.15 


4.1 


44.5 


7.24 


3.76 


0.31 


3.04 


0.70 


1.59 


0.36 


9.0 


3.74 


4.7 


48.6 


8.79 


3.76 


0.35 


3.05 


0.81 


1.58 


0.42 


10. 


7.07 


5.5 


45.9 


9.47 


3.65 


0.41 


3.00 


0.93 


1.60 


0.47 


11. 


5.46 


5.8 


44.3 


9.26 


3.59 


0.43 


3.05 


0.98 


1.61 


0.50 


12. 


3.96 


5.9 


48.2 


9.95 


3.37 


0.44 


3.09 


1.01 


1.65 


0.50 


13. 


4.87 


6.2 


51.5 


11.1 


3.22 


0.48 


3.04 


1.08 


1.65 


0.54 


14. 


6.54 


6.4 


48.2 


10.7 


3.15 


0.52 


2.91 


1.15 


1.69 


0.59 


15. 


7.71 


6.7 


48.3 


11.4 


3.07 


0.54 


2.82 


1.19 


1.77 


0.61 


16. 


5.39 


7.1 


48.0 


11.8 


3.08 


0.54 


2.86 


1.24 


1.84 


0.62 
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Table Bl. Continued. 



r 


V 


eV 


a 


eo" 


H/3 


eH/3 


Fe5015 


eFe5()15 


Mgb 


eMgb 


(") 


(km/s) 


(km/s) 


(km/s) 


(km/s) 


(A) 


(A) 


(A) 


(A) 


(A) 


(A) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


ID 0650 


0.0 


0.00 


1.9 


34.9 


3.39 


2.99 


0.20 


4.40 


0.40 


2.12 


0.24 


1.0 


6.43 


2.4 


39.4 


4.53 


2.43 


0.26 


4.16 


0.55 


2.08 


0.28 


2.0 


6.86 


2.9 


38.8 


5.38 


2.27 


0.31 


3.96 


0.65 


1.95 


0.33 


3.0 


9.83 


3.4 


39.9 


5.93 


2.14 


0.36 


3.58 


0.76 


1.97 


0.38 


4.0 


10.0 


3.8 


45.1 


7.20 


2.10 


0.41 


3.79 


0.89 


1.90 


0.45 


5.0 


11.0 


4.3 


45.4 


8.11 


2.04 


0.47 


3.68 


1.02 


1.94 


0.51 


6.0 


5.86 


5.0 


43.7 


8.24 


1.83 


0.52 


3.41 


1.11 


1.82 


0.55 


7.0 


6.88 


5.3 


41.5 


8.49 


1.80 


0.57 


3.14 


1.21 


1.83 


0.60 


8.0 


3.33 


6.0 


41.2 


8.59 


1.78 


0.63 


3.02 


1.33 


1.73 


0.66 


9.0 


2.06 


6.6 


40.6 


9.02 


1.78 


0.70 


2.84 


1.46 


1.62 


0.73 


10. 


5.57 


7.1 


40.3 


9.28 


1.76 


0.76 


2.73 


1.57 


1.55 


0.80 


ID 09 18 


0.0 


0.00 


0.6 


77.0 


1.02 


2.09 


0.07 


5.24 


0.17 


3.11 


0.08 


1.0 


19.0 


0.8 


77.3 


1.50 


2.05 


0.10 


4.87 


0.21 


2.97 


0.11 


2.0 


28.4 


1.2 


76.4 


2.21 


2.08 


0.14 


4.32 


0.30 


2.68 


0.14 


3.0 


23.4 


1.6 


75.7 


3.07 


2.07 


0.18 


4.05 


0.39 


2.54 


0.19 


4.0 


20.7 


2.1 


73.3 


4.11 


1.97 


0.23 


3.88 


0.47 


2.44 


0.24 


5.0 


17.8 


2.8 


72.3 


5.41 


2.01 


0.29 


3.63 


0.63 


2.30 


0.31 


6.0 


17.5 


3.7 


74.5 


6.92 


2.04 


0.36 


3.42 


0.78 


2.19 


0.38 


7.0 


9.35 


4.7 


75.9 


8.47 


1.93 


0.44 


3.44 


0.93 


2.00 


0.46 


8.0 


1.50 


5.8 


70.9 


9.87 


1.94 


0.53 


3.34 


1.13 


1.94 


0.57 


9.0 


0.99 


6.5 


72.8 


11.7 


1.97 


0.56 


3.24 


1.20 


1.86 


0.60 


10. 


6.39 


6.9 


72.2 


12.6 


1.91 


0.60 


3.26 


1.28 


1.82 


0.63 


11. 


9.30 


7.2 


70.7 


13.5 


1.87 


0.65 


3.43 


1.35 


1.77 


0.67 


12. 


8.84 


7.6 


71.1 


14.3 


1.78 


0.70 


3.23 


1.40 


1.80 


0.70 
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